Critical Lattice Size Limit for Synchronized Chaotic State in 1-D and 2-D Diffusively 

Coupled Map Lattices 
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We consider diffusively coupled map lattices with P neighbors (where P is arbitrary) and study 
the stability of synchronized state. We show that there exists a critical lattice size beyond which 
the synchronized state is unstable. This generalizes earlier results for nearest neighbor coupling. 
We confirm the analytical results by performing numerical simulations on coupled map lattices 
with logistic map at each node. The above analysis is also extended to 2-dimensional P-neighbor 
diffusively coupled map lattices. 
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I. INTRODUCTION 

In recent years, synchronization of coupled dynamical 
systems [l|, H H| has become an important area of re- 
search for their applications in a variety of fields includ- 
ing secure communications, cryptography, optics, neural 
networks, pattern formation, geophysics and population 
dynamics [H, HUB]- In particular, the stability of syn- 
chronized state in coupled map lattices (CML) with vari- 
ous coupling schemes have been studied extensively [E H 

EE EE EI El 0, EE EE EE EE El IMEiillll- Tobe 

specific the CML with diffusively coupling has attracted 
considerable attention in recent studies. In such systems, 
the synchronized state is not stable when the number of 
nodes exceeds a certain critical limit and each node is 
coupled only with its nearest neighbor [IE EH • In this 
paper, using the formalism put forth in Refs. [IE EH 
we derive an exact analytic expression for this limit for 
more general case of P-neighbor coupling. Further, the 
results are verified through numerical simulations in cou- 
pled logistic map lattices. All the analysis are carried 
out in both ID and 2D CMLs. Studies similar to our 
present work, but for coupled oscillators, are reported in 

Rcf. MM. 



II. CRITICAL SIZE LIMIT IN ID CASE 

Consider 1-dimensional coupled map lattices with P- 
neighbor diffusive coupling represented by 

I p 

xj-(n + l) = /(xj(n)) + ynXl a p [/( x j-p( n )) + 



p=i 



/(x,- +p (n))-2/(x i (n))], (1) 

where Xj is a Af-dimensional state vector, j represents 
the lattice site, L is the lattice size, a p is the coupling 
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strength between jth map and its pth neighbor, and the 
evolution of the map at jth site is described by /(xj (n)) . 
Also periodic boundary condition is imposed and the syn- 
chronized state (synchronization manifold) is defined by 
xi(n) = X2(n) = ••■ = Xi(n) = x(n). Since a p is a 
very general coupling coefficient, the long range model 
proposed by Antenedo Q can be incorporated into the 
above equation. 

Linearizing (fT]) around x and performing 
the discrete spatial Fourier transform i]i(n) — 
jr^2j_ 1 ex.p(—i2Trjl/L)zj(n), the resulting form af- 
ter simplification (see Refs. [E E3 f° r details) is 



= hi + In 



1 - — ^2a p sm 2 (TTpl/L) 
%= 1,2, ...,M; 1 = 0, -1. 



(2) 



Here /Xj(Z)'s are the Lyapunov exponents corresponding 
to Ith mode and h^s are the Lyapunov exponents of the 
isolated map ordered as hi > h% > ■ ■ ■ > hyi- The 
mode I = corresponds to the synchronized state and the 
other modes represents its transverse variations. Hence 
Hi (I) gives the largest transverse Lyapunov exponent for 
the mode I ^ 0. Therefore the stability of synchronized 
state is ensured if Hi{l) < for all I ^ 0. However, 
the symmetry in Fourier modes reduces this condition as 
< for I = l,2,...,i/2 ((£ - l)/2 if L is odd). 
Thus the stability condition reduces to 



2 

1 - — ^2a p sm 2 (irpl/L) 
p=l 

1 = 1,2, 



< exp(-/u), (3) 
L/2 or (L - l)/2, 



and this expression is also obtained in Ref. p. We use 
this condition to derive the expression for the critical 
lattice size limit in the rest of this paper. 

Let A; = 1 — -p Y,p=i Opsin 2 (npl / L) and define A max = 
max{A;}, A m j n = min{A/}. Then the above stability con- 
dition can be rewritten as: 
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Amax < exp(-/ii), A min > 



exp(-/ii) 



(4) 
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FIG. 1: (a) The variation of Li,p with hi, (b) the variation of 
L\,p with P, (c) the variation of e c with P in one dimensional 
coupled logistic map lattices, obtained for r — 1.9, and (d) the 
variation of Li,p with P in one dimensional coupled logistic 
map lattices, obtained for r = 1.9 (hi = 0.5554). 



Let Af = 1 — ^EpLi™ 2 ^?'/^)! e = miri{a p }, and 



as A* lax = max{Af}. 



define an upper bound on A n 
Therefore the first stability con 

cxp(-hi). Similarly, let \( = 1 - 2f J2p=i sm^^pl/L), 
e' = max{a p } and define a lower bound on A m i n as 
A* nin = min{A; }. Hence the second stability condition 
in |4]) is ensured if Aj^ in > — exp(— hi). 

We start with the simplest case of nearest neighbor 
coupling for which results already exist [Til . [l~3| and then 
we extend the idea to a more general P-neighbor cou- 
pling. In the nearest neighbor coupling case, Aj^ ax = 
1 — 2e sin 2 (7r/L). Therefore the first stability condition 
is satisfied if e > \ s ^F^~^q , where e < a p Vp. Conse- 
quently, in terms of the coupling coefficients, the above 
condition becomes 



a„ > 



1 — exp(— hi) 
2sin 2 (7r/i) 



Vp. 



(5) 



Similarly, A* 



1 — 2e. Hence the second stability con- 



dition is satisfied if e' < 



l+cxp(-/ti) 



where e' > a p Vp. In 



terms of coupling coefficients, we get 
1 + exp(-/ii) 



yp. 



(0) 



Combining the two inequalities ©fe©, we get the well- 
known [ill Il3j final stability condition as: 



1 — exp(— hi) 



< a„ < 



1 + exp(—hi) 



Vp. (7) 



2sin 2 (7r/L) " " 2 

As L becomes larger, the above stability range becomes 
smaller and at a particular critical value of L, the range 
shrinks to zero. Beyond this critical value of L, the stabil- 
ity condition J7]) is violated and hence the synchronized 
state can never be stable. At this critical value of L, one 
can replace the inequalities by equality signs in Eq.([7]) 
and get 



L 



i.i 



Int 



- 1 (y / tanh(/i 1 /2)) 



(8) 



where Li.i is the maximum lattice size that can support 
synchronized chaos in an one-dimensional nearest neigh- 
bor diffusively coupled map lattice. 

Let us now turn to the more general case of P- 
neighbors diffusively coupled map lattices for which no 
previous analytical results exist. However, similar results 
do exist for coupled oscillators [3, HH but they are in 
terms of numerically computed stability ranges whereas 
we give analytical expressions. After making use of some 
simple trigonometric relations, the expressions for A^ and 



Af take the forms 



Af = 1 



1 - 



and 



1 - 



sm(PTrl/L)cos((P + l)nl/L) 
Psin(vr//i) 



sm(P7rZ/£)cos((P + l)nl/L) 



(9) 



Psm(nl/L) 



(10) 



respectively, where e < a p < e' Vp. The expression inside 
the square bracket takes its lowest value when / = 1 and 
it takes its highest value for the mode / = lh =Int[Li i/2], 
for all values of P. Following the same procedure as in 
the nearest neighbor case, we finally get 



1 — exp(-hi) 

a(Pv/L)cos((P+l)ir/L) 



Psin(7r/L) 



1 + exp(-/ii) 

i(Pwl h /L)cos((P+l)TTl h /L) 
Psin(7Tl h /L) 



Vp. 



(11) 



At this critical coupling strength e c (a p =e c Vp) the ex- tremes values of a p coincide and the above expression 



3 



becomes 



sin(P7r/Li iP )cos((P + 1)-k/L 1>p ) 
sin(7r/Li ) p) 



P 



1 - tanh(/ii/2) 1 



where It = Int[Li i i/2] and L\ t \ is given in Eq. ©. The 
critical lattice size limit Li.p is obtained by solving the 
above transcendental equation numerically. In a special 
case of P = L/2 (or [L — l)/2 for odd L) we get /ii = 
2tanh~ 1 (l). This result indicates that the synchronized 
state is always possible for globally coupled map lattices 
as long as hi < oo. The dependence of L\.p on the 
maximum Lyapunov exponent of the isolated map (hi) 
and the number of neighbors coupled (P) are shown in 
Fig. QJa) andQJb). It is observed that Lip increases 
almost linearly with P for a particular value of hi, and 
decays with hi for a particular value of P. Also, all 



J 



sin(P7rZft/Lip)cos((P + 1)^/Lip) 
Psin(7r// l /Li ) p) 



(12) 



r 



the results are confirmed numerically by considering the 
logistic map (defined by x(n + 1) = 1 — r[x(n)] 2 ) at each 
node. The variations of the critical coupling strength e c 
and the critical size limit Lip with P (for r = 1.9) are 
shown in Fig. [He) andQJd). 



III. CRITICAL SIZE LIMIT IN 2D CASE 

Now we consider 2-dimensional coupled map lattices 
with P-neighbor diffusive coupling of the form 



Xj, fc 



1 p 

(n + 1) = /(xj, fc (n)) + — { a P [f(*]-pA n )) + f(xj+P,k{n)) - 2 /( x j,fc(«))] + 

[/(x j;fc _p(n)) +/(xj )fc+p (n)) -2/(x 3 -, fe (n))] J, 



p=i 



r 



(13) 



where x^fc is a M-dimensional state vector, (j, fc) repre- 
sents the lattice sites, L is the lattice size, and a p and b p 
are the coupling strengths between (j, k)th map and its 
pth neighbor along j and k directions, respectively. 

In this case, the stability condition for synchronized 
state is 



! p 

1 — — [apShi 2 (7rpZ/L) + 6pSin 2 (7rprn/Z/)] 



< exp(-/i x ), (14) 



where Z, m = 0, 1, . . . , L - 1, (l,m) ^ (0, 0). 

If we define !/; >m = 1 - p^ =1 [a p sin 2 (7rpZ/L) + 
& p sin 2 (7rpm/L)] , (Z,m) ^ (0,0), ^ max = max{i^, m }, and 



i^min = min{^ im } then the above stability condition be- 
comes: 

fmax < exp(-hi), f mi „ > -exp(-Zii). (15) 

Performing an analysis similar to the ID case, we ob- 
tain the expression for critical lattice size limit for nearest 
neighbor coupling as 



L 2 ,i = Int 



(16) 



sin _1 (^2tanh(/n/2)) 

In the case of P-neighbor diffusive coupling, we obtain 
the expression for the critical lattice size limit L 2 ,p as 



sin 



(Pn/L 2 .p)cos((P+l)TT/L 2 ,p) 



sin(7r/L 2 ,p) 



= P 



1 - 2tanh(/ii/2) 1 - 



i(Pnl h /L 2 ,p)cos((P + l)nl h /L 2i i 
Psm(irl h /L 2 .p) 



(17) 



where lh =Int [£2,1/2], £2.1 is the critical lattice size pled map lattices and is given in Eq. (|16[) . L 2 ,p is ob- 
limit for 2-dimensional nearest neighbor diffusively cou- 
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FIG. 2: (a) The variation of L2.P with hi, (b) the variation of 
1/2, p with P, (c) the variation of e c with P in 2- dimensional 
coupled logistic map lattices, obtained for r = 1.5, and (d) 
the variation of Li,p with P in 2-dimensional coupled logistic 
map lattices, obtained for r = 1.5 (hi = 0.2378). 



node. The variations of the critical coupling strength e c 
and the critical size limit L 2t p (for r — 1.5) with P are 
shown in Fig. [^c) and [2d). 



IV. CONCLUSIONS 

We have presented expressions for the critical lattice 
size limits (L% p and L 2j p) for both 1 and 2-dimensional 
coupled map lattices with P neighbor coupling. In both 
the cases, the value of these critical size limits increase 
almost linearly with the number of coupled neighbors. 
In addition, all the above results were verified through 
numerical studies using coupled logistic map lattices. 
Moreover, as P increases to the global coupling limit we 
showed explicitly that the critical size limit tends to infin- 
ity. However, our results are not valid for discontinuous 
maps (for example, tent and Bernoulli maps considered 
111 Refs. ®M) 

since the linear stability analysis used in 
this study requires the maps to be continuous. 
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